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SOME PROBLEMS 
REGARDING THE INFLUENCES 
OF THE MATERIAL CHARACTERISTICS 
ON THE TEMPERATURE DISTRIBUTION 
IN A DISC BRAKE 


Stefan VOLOACA,' Gheorghe FRATILA” 


Abstract. The paper continues previous researches of authors, where infrared images of 
brake disc were taken in real time experiments. A comparison between real results and 
the results obtained from finite element analyze in ABAQUS was made. It is observed that 
material properties of brake disc influence the thermal shock resistance. The 
thermograms analyze, obtained by experiments, makes possible to highlight the 
conductivity which has an essential influence over thermal stress. 
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1. Introduction 


Brake disc architecture design and optimization activities are followed by 
theoretical simulation and experiments based on a testing methodology with its 
accuracy [1]. 


As an example ABAQUS (as a software for finite element analyze), infrared 
thermal camera (as a precise experimental acquisition device of thermal 
distribution) and histograms of thermo grams (as a testing methodology). 


Every research activity, in simulation and experiments, has the aim to study the 
phenomenological and analytical of braking in order to obtain new information 
regarding the resulting thermal energy from disc brake (figure | and 2). 


For a high efficiency of brake disc, the resulted heat must be fast released, taking 
into account the cumulative rise of temperature [2], having unwanted effects over 
thermal stress thus resulted [3]. 


'University “Politehnica" of Bucharest, Automotive Engineering Department, Faculty of Transports 
(voloacastefan @ yahoo.com). 

"University “Politehnica" of Bucharest, Automotive Engineering Department, Faculty of Transports 
(ghe_fratila@ yahoo.com). 


126 Stefan Voloaca, Gheorghe Fratila 


104 Second braking 


a 


4 lia braking 

__ 

100 150 200 250 300 350 
Disc brake temperature [°C] 


Thermal stress [MPa] 
[o> 


Fig. 1. Intensive braking influence over disc Fig. 2. The phenomenological principle of 
brake temperature. Thermal stress variation thermal energy dissipation — radiation, 
after two brakings can be linear expressed, as convection and conduction 


results of paper [2]. 


It is mentioned, knowing thermal stress impact but having a small relevance than 
knowing their effects, design a durable disc brake. Forward is important to study 
the thermal shock as one of the significant effect of thermal stress. The thermal 
shock is a result of rapid temperature change during braking process, especially 
due to low thermal conductivity of a brake disc. The thermal stress resistance Rr 
can be written as [4]: 


~ ko,(1—-v) 
(ona 


R 


i (1) 
where k is thermal conductivity, or is maximal tension that material can resist, o - the 
thermal expansion coefficient , E - Young's modulus and v is the Poisson ratio. 


Observing relation (1), the thermal shock depends directly on conductivity. Cast 
iron is a common used metal for brake discs production. Its conductivity 
coefficient decreases with the rise of temperature [5, 6] having the effect of 
blocking the heat inside the disc further causing thermal shock. The quantitative 
evolution of thermal shock starts from heat balance corresponding to an 
elementary period of time dt of friction between brake disc and brake pads. Heat 
balance is expressed by relation: 


Q=0,+0,+0, (2) 


where: Q is the quantity of thermal energy as a result of braking, Q; — quantity of 
energy released by convection, Q2— quantity of energy released by radiation, Q3 — 
quantity of energy dissipated by conduction inside the brake disc. 


The risk of thermal shock is higher with a high heat quantity not dissipated. Due 
to the importance of the equation (2) further analyses of thermal properties of 
brake disc material are made: 
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e thermal conductivity k[W/(m-K)], defined as the amount of heat flux 


that would pass through a certain material depending on the temperature 
gradient over that material and can be written [6]: 


k=a-p-c, (3) 

where o — thermal diffusivity, p[kg/m*] - material density and c, — specific 
heat; 
e specific heat c, [J/(kg-K)], an intensive property which means that it is 
independent of the mass of a substance, defined as the amount of heat 
required to. raise the temperature of one gram of a substance by one Celsius 
degree; for cast iron, after [7] c = 420 J/(kg-K); 


e thermal diffusivity a [m7/s] is defined as the ratio of thermal conductivity 
to heat capacity; after [6] cast iron has a= (0.12...0.18)-10~“[m? /s]; 


© emissivity ¢ of a material is defined as the ratio of energy radiated to 
energy radiated by a black body at the same temperature. For polished cast 
iron, after [7], for a temperature of 24 °C, the emissivity is 6(24°C) = 0.64. 

A connection between first three characteristics is given by relation [2, 6]: 

— (4) 
pC, 
Forward, the authors will establish some relations between all four material 
characteristics. 


2. Influence of material properties over heat balance 


2.1. The influence of convection over the heat balance 


Quantity of thermal energy Q; released to environment by convection can be 
written [7]: 

Ol ng A. FE -T,) (5) 
where: 
meq - the medium convection heat transfer coefficient [W/(m? K)] into the 
[T...7,] and A, - the surface area of the disc brake [m7]. 


mD° _ 70.25)" 


As an example, for (T - T,) = (250-20) °C; A, = ji 


=(0.05 m’, fora 


variation of /,,. =10...100 W/(m’K), results 


Q, = (10...100) -0.05 - 230 =115...1150 W; 
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admitting 


=50 W/(m’-K), 


environment air 
results Q, =575 W the thermal energy transformed in time interval dz. 


Thermal energy released by convection is not dependent by material properties; it 
depends on dimension and geometry of brake disc. 


Facts that influence the heat release by radiation and dissipation by conduction — a 
bigger radius — a bigger mass m and a good dissipation by conduction: 


Q,=m-c, FT) =k: A, —-T,)- 8 (6) 
where g — brake disc thickness. 
Noyes and Vickers with their experiments [7, p. 27] concluded, at that time, the 
loss of heat by convection is significant (60%) 
2.2. Radiation influence over brake disc heat balance 


Thermal energy quantity Q»2, dissipated, depends on local temperature T and 
emissivity coefficient s. 


Overall, total thermal energy radiated Q> [Watts] in a time period df can be written 
due to Boltzmann’s equation: 


Q,=s-0-A-T* (7) 
where o = 5.67x10 * Wm -K “ is Stefan-Boltzmann constant; 
For a temperature 
mD* _ ™(0.25)" 
4 4 


results the thermal energy in time period dt 


T= 250°C; Af> = 0.05 m’* and ¢=0.4 


Q, =0.4-5.67-10 -0.05-250* =88 W. 


Hence, by radiation, a small part of thermal energy is dissipated, the rest of it is 
accumulated in disc and without conduction dissipation thermal stress appears. 


The material properties and the machinability of brake disc, have influence over 
the incertitude of intuitive estimated studying figures 3-6. 


Figure 5 presents images and their histograms for three values of emissivity 
1, 0.7, 0.5. Figure 6 presents the variation of luminosity with emissivity (for each 
histogram’s mean). 
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Fig. 3. The influence of radiated thermal energy of brake disc (expressed by its radiance) for two 
temperature values a) 100°C, b) 200 °C and four emissivity values (1; 0,5; 0,3 and 0,2). Both 
figures show the necessity of correct emissivity estimation, bigger errors of emitted radiance 
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appear at high values of brake disc temperatures (c). 
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Fig. 4. Correct emissivity setting and material state influence over thermo grams. With a high 
roughness, the temperature variation is smaller and relatively uniform. As a conclusion, at high 
roughness, thermal stress for different areas of brake disc is smaller. 
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Fig. 5. Images and their histograms with three set values of emissivity (1; 0.7; 0,5). 
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Fig. 6. The variation of luminosity with brake disc emissivity 


2.3. The influence of conduction over heat balance 


Conduction represents the heat transfer by means of molecular agitation within a 
material without any motion of the material as a whole. If one end of a disc brake 
has a high temperature, then energy will be transferred down to the hub which is 
colder. The first end has a high speed of particles that will collide with the slower 
ones with a net transfer of energy to the slower ones. 


Thermal energy Q3 dissipated by conduction is quantified by thermal conductivity 
coefficient k. The coefficient depends, too, on material characteristics and disc 
brake treatments. 
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For heat transfer between the two plane surfaces of the disc brake, such as heat 
loss through the wall of a house, the rate of conduction heat transfer is: 


Q, = KAT ot = Tia) (8) 
t & 


where A-area, d-thickness of barrier, Q3-heat transferred in time t. Total thermal 
energy Q3 in a time period can be expressed by relation: 


0 eae? (9) 
mD° — 70.25) 
\ TA 


An example is: T= 250°C; A, = =0.05 m*taking into account 
relation (2), results: 

k=a-p-c, =0.3-10°1/K-1800 kg/m? -420J/kg-K =0.22 
More over Q, = 0.22-0.05-420=88 W 


For a high temperature, the accumulated heat is released more by radiation and 
less by conduction, than a low temperature. Once with the temperature the 
molecular event grows too and the conductivity coefficient comes down. Thus, 
figure 7 presents a graphic resulted from experiment results in paper [6]. 
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Fig. 7. Thermal conductivity as a function of temperature for a cast iron, after [6]. 


For short period of braking (4+6 seconds) a mean value of conductivity can be 
considered (49...50). 


The dependence of k coefficient, on material characteristics and brake disc 
machinability, can be determined from diagrams of figure 8. This figure presents 
temperature variation in time for two types of cooling for vented brake disc. This 
temperature variation was made to determine emissivity variation after a 
methodology presented by him in paper [8]. Pedrag and al. in paper [1] used a 
similar methodology to determine conduction coefficient. 
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Fig. 8. Temperature variation for two types of cooling. 


For cast iron , thermal conductivity k is quite high, and between usual temperature 
met at first braking (under 250°C) varies between 40...50, as a function of 
temperature value 7. At a given temperature, thermal conductivity is proportional to 
the raise of temperature. This behavior is quantified in the Wiedemann-Franz Law: 

k=o,°L-¥ (10) 
where Lorentz’s number L is 2.44x10° WQK?® and o¢ = electric conductivity of 
brake disc material. 


3. Approach the ideas 


3.1. Supplementary ventilation of brake disc 
Thermal stress is defined by relation: 
o=E-a-AT (11) 
It’s more important its time variation: 
2 =E ee -AT+a@ 


Due to small braking period (approximately 4 seconds) high thermal conductivity 
and emissivity modifications occurs. 


“ ) (12) 


From figure 7, taking into account relation [9]: 


k=“ [int, -Int,]/[Z,-T] 
4n 
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It is established 
q n20,5+-In5,5 9 q 4.7 


k=“ fine, —int/(0, -T) = a 500 
Ar ° 4x 101,9-101 Ar 0.9 Ar 


For high ventilation (forced) results: 


qd q In20,5+1In5,5 q 4.7 qd 
K ventilation = ——Un tp —n 4 A cd = = 0,36-— 
ventilation = 47m — A Td 74-101 40 13 An 
The ratio 
& ALA _ lg. 
Rventiltion 0.36 


A strong decrease of conductivity appears in case of forced cooling by ventilation. 
The temperature rise can be a phenomenological explication of this, for the same 
period of time. 


A rise of conductivity was expected but the molecular event was high due to high 
rise of temperature and the effect is reverse — a lower conductivity. 


For the vented brake disc, the conduction is smaller, local temperature is higher so 
the heat will be dissipated by radiation. 

woud -7* 3 
Of ATE. 4, =-o(=const)- ail 
Q, k-A,-T 


An analyze of the ratio variation was made: 


From equations (6), (7) and (9): 


e For a thick disc (a bigger mass) the coefficient k is smaller for the same 
quantity of thermal energy Q3, see equation (3); 


e For a high temperature, the coefficient ¢ is smaller for the same 
accumulated thermal energy Q; see equation (7); 


const. 


e From figure 7. and due to ratio sa TIK} for T €[50...250],k ¢[52...47], 


the coefficient k = 53-0.025 T; 


4 3 
21s = .o(= const): = us 
Q,  k- AT 
braking from 100 km/h ... 50 km/h (where brake disc temperature is under 


250°C), the emissivity and conduction varies in opposite ways by 
temperature (figure 9). 


, for intensive 


e For constant ratios 
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Fig. 9. The ratio variation 1 in relative units reported to temperature variation. 


Similar variation of conductivity with temperature is given for steel [10, 11]. 


Must be mentioned, from paper [12] the high effect of thermal stress appears in 
the first 4 seconds of the braking. 


3.2. Numerical and experimental analyze. Results interpretation 


Nowadays traffic is characterized by a high density of motor vehicles and 
intensive braking, especially at town entrances, from usual speeds of 100 km/h ... 
60 km/h. 


Due this, a disc brake must be efficient and must permit a fast stop of motor 
vehicle in 2+6 seconds [2] from a speed of 100 km/h. An example is given by a 
braking period of 5.48 seconds for a Volvo motor vehicle at a temperature of 
27°C [2]. Considering the high effect of thermal stress in the first seconds of 
braking [12], the authors in their finite element analyze considered a braking 
period of 4 seconds. 


Due to previous results of experiments made and previously published by the 
authors, the research continued with a finite element analyze, made in ABAQUS. 


The aim of these was to validate an interrelation between emissivity experimental 
determined and thermal simulation of conductivity with mentioned software. 


To validate results, a comparison between images obtained from experiments and 
simulation was made. Thermal images were processed with MikroView [13] and 
Termografie [14]. The comparison was made for two velocity values 80 km/h and 
50 km/h (figure 9, 10 and 11). Temperature repartition is non-uniform in radial 
direction due to high convection and radiation to surface extremities. 
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a. analyse with ABAQUS Vit b. analyse with MikroView 
Fig. 9. Temperature repartition over disc surface at 80 km/h. 


a. analyse with ABAQUS of temperature in elements b. analyze with Termografie 
— the selected element temperature is 143 °C 


Fig. 11. Temperature repartition over disc surface at 50 km/h. 


Figure 11.a presents an element of modeled brake disc which has a temperature of 
143°C. The thermal image (Fig. 11.b) was processed with software named 
Termografie which made possible to select an area identical as position to one 
obtained by simulation. For that area the temperature was 147 °C. 


The results of this comparison are in good concordance (the error between 


experiment and simulation results are about ¢ = ali 100 = 2,72%). 
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Fig. 12. Thermal energy transfer due to radiation and conduction as a function of thickness. 


4. Conclusions 


1. 


Thermal properties of brake disc material are important in establishing 
braking efficiency. 


Convection influences thermal stress by brake disc radius - the bigger it is, 
the limit thermal layer is higher and thermal bridging takes place at the 
disc periphery. 


Conduction is significant influenced by temperature evolution for an 
intensive braking; the conduction coefficient comes down with 
temperature and emissivity rise. 


Having emissivity variation is possible to know conductivity variation in 
heat balance due to dependence between appreciation errors of both 
characteristics; thus the conductivity variation can be followed like the 
emissivity variation. 


Emissivity variation at a moment can be highlighted by luminosity of 
thermal image; the same variation of luminosity permits the evaluation of 
brake disc material’s conductivity and forward to thermal stress. 
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